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PARTICLE BED REACTOR 
MODELING 


■ PRESENT THERMAL- HYDRAULIC 
SYSTEM MODELING TOOLS B&W USES 
FOR NTP SYSTEMS 

■ FOCUS ON PARTICLE BED REACTOR 
TECHNOLOGY AND THERMAL 
HYDRAULIC METHODS. 


View Graph 4 - Particle Bed Reactor Modeling 


1 . ^ttonwl-hydnulic .yfem *»b UW urn fcr nucta, th*™. 

^ed reactor technology and the thermal-hydraulic methods used I 
*n who feel that thermal .hvHraiilu* nwiai;.. i • 
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NTP: Systems Modeling 


582 


NP-TIM-92 



PBR CORE CROSS SECTION 



View Graph 5 - Radial Cross Section of Particle Bod Reactor 

ThU view graph shows a radial cross section view of the reactor system we will be discussing today. This system is a generic 
particle bed reactor syliem made up of 37 fuel elements as shown by the red circles. The blue area surrounding the fuel elements 
are hexagonal moderator blocks. Some of the holes shown in the blocks are for propellant flow through the moderator. 

This core is surrounded by a reflector and twelve control drums which are in turn surrounded by a pressure vessel. Details 
of particle bed reactor systems were presented in several papers at this workshop and won't be covered here. 
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PBR BLEED CYCLE 



View Graph 6 - PBR Bleed Cycle 

~ -Sr !*££?£>%;££%%£ l r r^“ *- 

y pknumi end have a single pan cooling loop. Depending on design requirements flow split and untie nut cmratf « 
on be used for any combination of the moderator, reflector or nozzle wall flow paths. * ooncqtU 

The propellant exit* the moderator and ii collected in a plenum above the core. It is then sent through the fuel mh 

mT'ZZ?' 1 *” 0 ' 1 * T«rt« outlet tempentures are nominally very high to maintain high KR^Mach nunAeTiTabout 

. ** pur P° sc * <* TT * c * or ""Wi ni there are three area* which are usually discussed separately since they retrain different *v«*« 

ta ** eMirt par ” ch> bed •»“*» "**« ty£*!. 

mmbUei. TOl qutm modeling will not be dnettued here today. The other two areas are fluid flow In the entrance and eak 
plenums of the reactor system and finally modeling of fluid flow through the particle bed fuel element. 
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FUEL AND MODERATOR FLOW 

PATHS 



View Graph 7 - Fuel Element Flow paths 

TWtisa view of apartide bed fuel dement with flow paths shown by arrows. The red hatched (outer) am ii the moderator 
•ection; the onnge vea it the Aid bad and the green treat the burner (hot) and outer (cold) frill that hold the fuel particles. 

i 1 ** 1 *“ “ *® moderl,or 10 cool it, then to a plenum at the entrance aide of the fuel element, or 

directly into the fuel element. Orifice* of the element can be done at either die moderator entrance or the fuel dement entrance 

The gas enters the oold frit which is at the outer annulus of the fuel element, then passes through the fuel bed and hot frit 
where it turns and flows out the outlet channel. ^ me me ueo, ano not ant 


Target outlet temperatures are high to maintain high specific tmphne. Mach number it appnniraaidy 0.25 at the outlet. 
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PBR MODELING REQUIREMENTS 


■ 1. FLUID FLOW THROUGH A PARTICLE 
BED 

■2. COMPRESSIBLE AND 
INCOMPRESSIBLE FLOW 
■3. SINGLE and TWO-PHASE FLOW 

■ 4. COUPLES FLUID FLOW and SOLID 
HEAT TRANSFER 


View Graph « - PBR Modelin» RequiraraenU 

flow Setos e. (fl-M tort, irafcpeodee.orpraourad^e.) S 

compfBttible flow or with i Kpetate treatment unng equations for thermally expandrtle flow. tt * U “ 0W wed for 
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PBR MODELING REQUIREMENTS, 

Cont. 


■ RANGE OF SINGLE TO 
MULTI-DIMENSIONAL MODELING 

■ TRANSIENT AND STEADY-STATE 
ANALYSIS 


View Graph 9 - PBR Modeling Requirements (Continued) 

The computer codes used to analyze the fuel element will need mulli -dimensional capabilities. The systems level analysis will 
use primarily one-dimensional techniques. Both transient and steady state analysis will be required to cover the wide range of 
operating and accident modes. 
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CHARACTERISTICS OF PRR and 
NTT MODELING 


■ 1. FUEL ELEMENT FLOW-TO-POWER 
MATCH 

■2. REACTOR FLOW-TO-POWER MATCH 

■ 3. BED TO COLD FRIT HEATING 
EFFECTS 


View Graph 10 - Characteristic! Of PBR and NTP Moddini 
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CHALLENGES FOR PBR and NTP 

MODELING 


■ 1. START UP TRANSIENTS 
-2. DECAY HEAT 

- 3. THROTTLING CONDITIONS 
■4. ACCIDENT TRANSIENTS 

■ 5. PRE-TEST PREDICTIONS 
■6. COMPONENT HEATING 


View Graph 11 - Challenges for PBR end NTP Modeling 

P M* yi HMber of applications of modeling required for i PBX reector. There eln include use of modeling 
fa drel gdn^wra md performing poet-test evaluations. Examples of system analyses for Decay Heel cooling end Shut Up Transient 
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THERMAL HYDRAULIC 
COMPUTER CODES 


■ 1. OTV ENGINE - B&W 

- PARTICLE BED FUEL ELEMENT DESIGN 
SPECIFIC 

■ 2. TEMPEST- BATTELLE NORTHWEST 

- GENERAL 3-D CFD ANALYSIS 
- 3. SAFSIM - SANDIA 

- NETWORK SYSTEMS ANALYSIS CODE 

■ 4. SINDA/SINFLO-NASA 

- DETAIL THERMAL ANALYZER 


View Onphs 12-21 - Thermal Hydraulic Computer Codes, Code Ctpebilidu end Limitation. 

vie " * how ** m hi OT thermal-hydraulic codes which have been used by MW for analysis of NTP 
^ wtdtwmw of dwir mpebiliti* arai llmltmlo*. -now do*.'. allow . WlTcZL ofd** ^ ^ 

one called OTV B^b*. This computer The pnmmy oode la lUa ctus w* 

coda are a good iadiartion of why a J»ree number of codes areuaed In *“ "^a**!*** 

n. „x* u* i43K«iSr^2nS2 
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CAPABILITIES FOR 
PBR/REACTOR APPLICATION 


■ OTV-ENGINE 

- PROVIDES "NOMINAL" FUEL ELEMENT 
DESIGN CONDITION 

- SPATIAL FUEL TEMPERATURE 
-PROVIDES "OFF-NOMINAL" DESIGN 
CONDITIONS 


THERMAL/HYDRAULIC CODES, 

COllt. 

■ 5. ANSYS - SWANSON, INC. 

- DETAIL THERMAL CODE FOR 
COMPONENT AND LOOP ANALYSIS 

-6. NEST-B&W 

- TRANSIENT ANALYSIS OF COUPLED 
NEUTRONICS, 

THERMAL-HYDRAULICS 

■7. ATHENA -INEL 

- 1-D TRANSIENT OR STEADY STATE 
SIMULATION OF SPACE REACTORS 
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CAPABILITIES, cont 


■TEMPEST 

- MULTI DIMENSIONAL CFD ANALYSIS 
-ALLOWS ANALYSIS OF ACTUAL DESIGN 
-ADDRESSES COMPLEX THERMAL/FLOW 
■SAFSIM 

-REACTOR AND ENGINE SYSTEM 
■ SINDA 

-GENERALIZED CONDUCTION AND 1-D 
CIRCUIT FIOW SPLIT MODELING 
CAPABILITY 


CAPABILITIES, Cont. 


■ANSYS 

- PERFORMS GENERALIZED DETAIL 
HEAT TRANSFER ANALYSIS 

~ PROVIDES GENERAL COUPLED 
FLOW/CONDUCTION HEAT TRANSFER 
FOR SPECIFIED (KNOWN) FLOW 
REGIONS 
■ NEST 

- EVALUATION OF SYSTEM CONTROL 
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LIMITATIONS 


■ OTV-E 

-STEADY STATE 

- NO REACTOR PHYSICS 

- NO CONDUCTION (gas or solid) 

- NO GENERAL FEATURE CAPABILITY 

- CHANNEL APPROACH TO FLOW (1-D) 


LIMITATIONS, cont. 


■TEMPEST 

- NO REACTOR PHYSICS 

- LIMITED TO ORTHOGONAL CURVELINEAR 
GEOMETRICS AT PRESENT 

- TIME STEP LIMITED TO "MATERIAL- 
COURANT 

■ SAFSIM 

-TIME STEP UMITED TO "MATERIAL- 
COURANT" 

- PSEUDO MULTIDIMENSIONAL (1-D FLOW, 
NETWORK HEAT TRANSFER) 
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LIMITATIONS, cont 


■ NEST 


-POINT KINETICS 
- QUASI -STEADY FLUID FLOW 
■ SINDA 


- MODEL DEFINITION IS TEDIOUS 

- FLOW IS INCOMPRESSIBLE 

- NO SPECIFIC PROVISION FOR FLUID FLOW 
THROUGH PARTICLE BED 

-STEADY STATE 


LIMITATIONS, cont. 


■ANSYS 

-STEADY STATE FLOW 

- INCOMPRESSIBLE FLOW ONLY 

- LACKS SPECIALIZED CORRELATION 
CAPABILITY (FRICTION, FILM 
COEFFICIENT, etc.) 

- PSEUDO MULTI-DIMENSIONAL (1-D 
FLOW, 3-D HEAT TRANSFER) 

■ALL CODES LISTED ARE SINGLE PHASE - 

WILL NEED TWO PHASE CAPABILITY 
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PHYSICAL CORRELATIONS 


■ SPECIFIC CORRELATIONS FOR 
PARTICLE BED 

- FILM COEFFICIENTS - ACHENBACH 

- FRICTION COEFFICIENT - ERGUN 

■ FUEL ELEMENT COMPONENTS (COLD 
& HOT FRITS) 

- MODIFY GENERALIZED 
CORRELATIONA FOR SPECIFIC 
APPLICATION BASED ON 
EXPERIMENTAL DATA 


View Graph 22 * Physical Correlations 

. . 7 *?. nC * l , t ? 0 v ! cw , provide some information on the second major component of systems modeling - the validity and 

determinallon of the physical parameters and condatkmi used for modeling of the system. These view graphs show well known 
oorretatkm^ have been used in particle bed modeling. They also identify the need for experimental verification of this data. 
B*w n*s performed many of the experiments required to verify this date. 
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MODIFIED CORRELATIONS 


■EXAMPLES 

- MODIFY ERGUN CORRELATION FOR 
COLD FRIT 

“ FRICTION FACTORS FOR BLOWING 
AND SUCTION FLOW 

- PARTICLE BED CONDUCTIVITY - 
ZEHNER AND BAUER 


View Graph 23 - Modified Correlations 

Con ? iklk ** th* h«vc b««i modified are shown here. Ttey include friction coefficients for cMd f riti. 
friction heton for Mowing and auction flow in dm entrance and exit annulus of the fod element and f y flrhr bed conductivity. 


NTP: Sy» terns Modeling 


596 


NP-TTM-92 



Comparison of Predicted Friction 
Factor and Exper imental Data 


1.6E+6 



4E+5 1 i * ■ . j_ I 

3E+5 5E+5 1E+6 2E+6 

Re/Dp (1/m) 


View Graph 24 - Comparison of Predicted Friction Factor And Experimental Data 

Thii view shows a comparison of a predicted friction factor correlation of a outer (cold) frit as compared to the design 
correlation determined from experimental data taken at BAW's Alliance Research Center, In this case, air was flowed through typical 
manufactured fritsand pressure drop measurements performed. This plot is a measure of the normalized friction factor as a function 
or Reynolds number. As you can see the design correlation, which has an accuracy of plus or minus 10%, is approximately 30 to 
40% higher than the theoretics] friction factor and shows a steeper increase with lower Reynolds number. 

. 1*1 addition *° ***** c °k* B&W has used experimental data for friction factors covering blowing and suction flow in 
the fuel element annulus and have plans for performing tests on particle bed conductivity. As shown on the previous view graph 
BiW currently uses the correlation of Zehner and Bauer for particle bed conductivity. This correlation was not developed for PBR 
applications and therefore will be experimentally verified. 
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FRIT PRESSURE DROP TESTING 
WITH H 2 , AIR, and N, 


T EST CONDITIONS 


p 

3.2 MPa 

T 

294 K 

(%) 

Re/Dp 


Sis 

tULizmi 

Kexn fio 

Air 

5.08 

5.51 

Air 

5.02 

5.36 

H 

5.38 

5.25 

H 

5.38 

5.27 

S 

5.04 

5.39 


Kexp-Kcalc 
KcaiCtlQ ) Kcalc 


5.43 

•1-1.5 

5.48 

-2.2 

5.19 

+1.2 

5.19 

+1.5 

5.47 

-1.5 


View Graph 25 . Prit Procure Drop Toting 
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COLD FRIT MASK FACTOR 

With Azimuthal Power Variations 



View Graph 27 - Cold Frit Muk Factor • Aiimullal Variation 

—t. 10 | — " ! n ooihe outer (cold) frit, fricdo. 

wen that da flow matcliat th* pm* to u*una| dw too ***„ iBnowoccwutetai^ 


an uaed lo design the frit 
qxtfs. 

The three curves tie for the 


k This view graph shows typical „ 

sarjssssRcsL'ir- 


comliU^inTp^W^* 1 — J* Prt— far d«r <«. (Mfa, oadUw) 


and start up 
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Decay Heat Flow Rate 


Average Flow Rate During Decay Heat Cooling 


Average Flow Rate During Decay Heat Cooling 




View Graph 28 - Decay Heal Flow Rale 

Thl* view graph shows a typical evaluation of the propellant flow rate required after shut-down to cool a panicle bed reactor 
under decay heating caused by the gamma and beta radiation being emitted by nuclear fuel after shut-down. The scenario used for 
decay heat cooling is to maintain a throttled overflow of propellant for approximately the first 100 seconds after shut-down to insure 
a cool geometry. The flow is gradually decreased to match the declining power output of the core until the 10% flow plateau is 
reached. This flow is maintained constant for a while due to stability considerations which I will discuss later. The system then 
converts to pulse cooling similar to that planned for the NBRVA engine. Pulse cooling continues through approxi matel y 10,000 
seconds or until the system gets to approximately one to two percent of full power. At this point a long-term closed cycle coding 
system would be used to keep the reactor cooled through some type of closed loop system. This system would radiate the small 
excess heat to space. The view graph on the right is a plot of the actual predicted flow to the optimum flow needed for this process. 
In this case optimum flow would be that flow needed to exactly match flow to system heat rale. As you can see there is a spike 
where the actual flow exceeds the optimum flow by approximately five times for a short period of time to accommodate instability 
limits. 


It should be noted that the numbers shown here were obtained with analysis of a single fuel element They do not account 
for flow splits in the total system. Also no mechanical analysis were performed to determine the effects of thermal cycling during 
pulsed cooling. 
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Stability Regimes 



View Qnph 29 • Stability Regimes 

This view graph shows two presentations of the same data. The one on the left using a log scale for the V axil ml the one 
Mthc i^ using a linear scale for the "x" axis. The m y m axis is a plot of a instability index developed by Bustard based on NBRVA 
*** and applied by Maim of Brookheven National Laboralory to particie bed systems. This index is the dfflermoe between (he inks 
and (he outlet temperature divided by the inlet temperature and here it is shown as a function of power density. If yoo focus vour 
attention to the view graph on the right, the opes area is that region where flow instabilities would not occur. The shaded ration 
ii where there are potential flow instabilities, In the case discussed here, the shaded area is only approached during decay heat 
cooling and Is not a factor in the operating regions. V ' 

I 11 !* on the left shows an example of how the unstable region shrinks as one performs more detailed analysis of 
flow instabilities. The curve shown at the right labeled 1-D stability limit is the analytical result <+**i~* ^ 9umu a 
cireka re present ■ ihift or the ooe-dimen»onal stability regime when inalyied with the computational fluid dyrmnks computer code 
TBMPEST. The darkly shaded areas show even farther movement when a particle bed system is aoalyred with tiueedlmntioml 
codre. In this case the area is shown shaded became no sharp boundary exists. Insteed we are predicting a gradually increestitt 
probability of flow maldistribution. The actual region of instability would have to be verified by experiment beca me of three 
uncertainties. These curves show the advantages of using multi-dimenstoaal analysts on these complex geometries. 


We need to note that this is not only a PBR problem - all gas reactors will need to accommodate 
flow/high delta T conditions. 


instability Limits at low 
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STARTUP TRANSIENT 
SIMULATION 



TIME(SEC) TIME(SEC) 


View Graph 30 - Start-up Transient Simulation 

Thi* view graph gives a representative example of an analysis performed for the start up of a particle bed reactor. This 
analysis was done with B&W*i NBST computer code system. It was performed lo evaluate the unusually high reactivity insertion 
{rofC [ flowing cold hydrogen during start-up of the system. In particular It was being used to evaluate the effectiveness of the control 
mechanisms to mitigate the large insertion of positive reactivity into the system during start-up. These slides show the percent power, 
percent hydrogen flow, hydrogen worth, and reactivity change of the system versus time over a period of approximately twelve 
seconds. This analysis shows the system can achieve and maintain design power. 

TTie start up scenario used here is "dry*. The reactor is taken critical before hydrogen flow is initiated. As hydrogen starts 
to flow one set of control elements Is moved to overcome the positive reactivity insertion by hydrogeo flow. Another set of 
control elements, with different characteristics from the first, is used to control power. The control algorithm controls to a demand 
Startup period white coostnM ty maximum power versus fiow requirements which are shown in this viewgraph. 
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PHILOSOPHY OF SYSTEMS 
MODELING 


■ THE PROOF OF THE PUDDING IS IN THE 
TESTING 

■ LEARN FROM EXPERIENCE 
- SKYLAB and HUBBLE 

■ SYSTEMS MODELING IS A GUIDE FOR 
PERFORMANCE AND TESTING. IT IS NOT 
THE FINAL WORD 


View Graph 31 - Philosophy of Systems Modeling 


-.rf,, ,, 7 ]” " » P ** erl l * tIltude ” fttomfkf *»wnl. ■wMIk tfati taunt U wW a wriN ^ 

mbV * m ****** «*< «*«««. <»». effa* tttt of phytfol pmM, 2d 


T * ie HuW> ** tetetcope had significant problems because it wasn't tested before launch. 
?!T other veJ,icici WM »!*»«*• This is not intended to pick on NASA, there 
taiea to tea. There were picked because they ere recent or more easily identified by NASA. 


Skylab was damaged during launch 
are other iad retries that hare similar 
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SUMMARY 


■ CHALLENGES OF PBR MODELING AND 
SYSTEM ANALYSIS 

■ COMPUTER CODES 

■ PHYSICAL CORRELATIONS 
■RESULTS OF ANALYSIS FOR DECAY HEAT 

COOLING AND STARTUP 

■ PHILOSOPHY OF SYSTEMS MODELING 


View Graph 32 - Summary 

In summary this presentation has covered the characteristics and some challenges of Particle Bed Reactor modeling. It covered 
the mqjor components of modeling; Computer codes, physical correlations used, a test philosophy, and selected results of decay heat 
cooling and start-up analyses. 

Finally, there was an appeal to sit of us to keep in mind the necessity of obtaining experimental data to verify systems 
performance and systems models. 
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FINAL THOUGHTS 


■NOBODY BELIEVES THE ANALYSIS 
EXCEPT THE ANALYST 

■ EVERYBODY BELIEVES THE 
EXPERIMENT EXCEPT THE 
EXPERIMENTALIST 

- Seen on NASA wall 

■ "PAPER REACTORS, REAL REACTORS" 
Admiral Hyman Rickover - 1953 


View Graph 33 - Pina] Thought* 

Sp«* * *” * eal ” * will po&zr ditrin| I rtcenl vuh Id iJie Hunttvilfc 

sxsstttsas: Era r. 
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Paper 

Reactors, 

Real 

Reactors 


Admiral Hyman Rickover, 
The Journal of Reactor 
Science and Engineering, 
June 1953 


An academic reactor or reactor plant almost always has the 
following basic characteristics: 1) It is simple. 2) Ms small. 3) It 
is cheap. 4) It Is light. 5) It can be bum very quickly. 6) It Is very 
flexible In purpose. 7) Very little development is required. It win 
use mostly off-the-shelf components. 8) The reactor is in the 
study phase. It is not being built now. 

On the other hand, a practical reactor plant can be 
distinguished by the following characteristics: 1) It is being built 
now. 2) It Is behind schedule. 3) It Is requiring an immense 
amount of development on apparently trivial items. Corrosion, In 
particular, is a problem. 4) It Is very expensive. 5) It takes a 
long time to bund because of the engineering development 
problems. 6) It Is large. 7) Ills heavy. 8) It is complicated. 

The tools of the academic reactor-designer are a piece of 
paper and a pencil with an eraser. If a mistake is made, it can 
always be erased and changed. If the practical-reactor designer 
errs, he wears the mistake around his neck; it cannot be 
erased. Everyone can see it. 


The academic-reactor designer is a dilettante. He has not 
had to assume any real responsibility in connection with his 
projects. He is free to luxuriate In elegant ideas, the practical 
shortcomings of which can be relegated to the category of 
“mere technical details.* The practical-reactor designer must 
live with these same technical detals. Although recalcitrant and 
awkward, they must be solved and cannot be put off until 
tommorrow. Their solution requires manpower, time and money. 

Unfortunately for those who must make far-reaching 
decisions without the benefit of an intimate knowledge of reactor 
technology, and unfortunately for the interested public, it is 
much easier to get the academic side of an issue than the 
practical side. For a large part those involved with the academic 
reactors have more inclination and time to present their ideas in 
reports and orally to those who will listen. Since they are 
innocently unaware of the real but hidden difficulties of their 
plans, they speak with great facility and confidence. Those 
involved with practical reactors, humbled by their experiences, 
speak less and worry more. 


Yet it is incumbent on those in high places to make wise 
decisions and it is reasonable and important that the public be 
correctly informed. It is consequently incumbent on all of us to 
state the facts as forthrightly as possible. 
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